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(Ul^he  steady  state,  zonally  averaged  circulation  of  the  middle  atmosphere  (15-125  km|  is 
studied  with  a quaaigeostrophic,  numerical  model  that  explicitly  includes  a self-consistent  calculation 
of  solar  radiative  heating  due  to  i6^nd^0^  absorption,  Newtonian  cooling,  Rayleigh  friction,  tropo- 
pauae  boundary  conditions  baseoW  cllmafological  averages,  and  the  effects  of  vertically  propagating 
planetary  waves.  We  find  the  direct,  radiatively  driven  pole  to  pole  circulation  at  solstice  is  sufficient 
to  account  for  the  cold  summer  menopause  and  warm  isothermal  winter  mesosphere  with  associated  - 
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20.  Abstract  iContinuedi 

tonal  )«ts  of  realistic  matnitude.  The  climatolof  ical  heat  and  momentum  nitxct  aaaociated  with 
planetar>'  wave  number  2 have  a neglifible  effect  on  the  mean  circulation.  With  planetar>'  wave 
number  1 no  steady  state  solution  could  be  obtained  due  to  the  formation  of  eaiterlles  and  hence 
critical  layers  in  the  winter  meeosphere.  We  also  find  that  the  radiative  haatinc  aaaociatad  with 
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THE  ZONALLY  AVERAGED  CIRCULATION 
OF  THE  MIDDLE  ATMOSPHERE 

I.  INTRODUCTION 

In  recent  years  there  has  been  considerable  interest  in  stratospheric  chemistry  as  a result 
of  concern  about  pollution  effects  on  the  ozone  layer.  In  the  lower  and  middle  stratosphere 
transport  processes  affect  the  ozone  density  distribution  as  a consequence  of  its  long  chemical 
lifetime.  The  motion  field  of  the  middle  atmosphere  is  thus  not  only  of  intrinsic  interest  but 
plays  an  important  role  in  our  chemical  understanding  of  this  region. 

Vertically  propagating  planetary  waves  forced  in  the  troposphere  provide  an  important  dynamical 
mechanism  to  couple  the  stratosphere  and  troposphere.  The  reflection  and  transmission  of  these  ultra- 
long  waves  are  strongly  dependent  on  the  structure  of  zonal  wind  jets  centered  near  the  stratopause 
(Schoeberl  and  Geller,  1977).  Further  progress  in  understanding  such  fundamental  dynamical  processes 
in  the  stratosphere  as  sudden  warmings  and  their  effects  on  the  lower  atmosphere  require  a thorough 
knowledge  of  the  dynamical  processes  that  control  the  polar  night  jet. 

The  dynamics  of  the  middle  atmosphere  is  often  described  in  terms  of  a zonally  averaged 
circulation  and  perturbations  from  this  average.  The  latter  include  a variety  of  wave  phenome- 
na which  are  principally  cont'olled  by  linear  processes  and  thus  amenable  to  accurate  quanta- 
tive  simulations.  The  symmetric  (zonally  averaged)  circulation  depends  on  the  nonlinear  eddy 
statistics  of  the  middle  atmosphere  through  the  divergences  of  the  heat  and  momentum  fluxes 
associated  with  the  wave  perturbations  and  on  solar  radiative  heating.  The  lack  of  information 
on  eddy  statistics  has  hampered  progress  since  the  initial  attempt  of  Leovy  (1964)  to  model  the 
ManuKiipt  lubmitted  July  13.  1977. 


symmetric  circulation  of  the  middle  atmosphere  At  solstice  Leovy  showed  that  radiative 


forcing  by  o2one  absorption  of  sunlight  would  drive  a pole  to  pole  mean  circulation  that  could 
explain  the  polar  night  jet,  the  cold  summer  mesosphere  and  the  warm,  isothermal,  winter 
mesosphere.  These  are  the  essential  observed  features  of  the  middle  atmosphere.  Since 
Leovy's  work  considerably  more  data  have  become  available  for  comparison  with  theoretical 
models.  Also  important  advances  have  been  made  m the  parameterization  of  infrared  radiative 
cooling  (Dickinson,  1973a).  accurate  representation  of  the  differential  radiative  heating  due  to 
insolation  (Park  and  London,  1974).  and  photochemical  acceleration  of  thermal  relaxation 
(Blake  and  Lindzen.  1973). 

In  this  paper  we  show  that  it  is  possible  to  numerically  simulate  certain  essential  features 
of  the  symmetric  circulation  with  accurate  input  parameters  and  representation  of  the  eddy 
heat  divergence  based  on  the  above  mentioned  advances.  While  our  approach  and  model  is 
similar  to  Leovy's  (1964).  we  solve  the  resultant  equations  numerically  and  relax  a number  of 
constraints  that  Leovy  imposed  to  produce  a mathematically  tractable  problem  for  analytic 
techniques.  We  allow  Newtonian  cooling  to  vary  with  height  and  latitude  whereas  our  only  ad- 
justable parameter.  Rayleigh  friction,  varies  with  height  only  for  lack  of  information  on  its  lati- 
tudinal variation.  In  addition,  we  have  attempted  to  include  the  contributions  of  planetary 
waves  with  zonal  wave  numbers  1 and  2 to  the  local  divergences  of  the  eddy  heat  and  momen- 
tum fluxes  by  explicitly  computing  their  wave  amplitude  with  the  model  of  Schoeberl  and  Gell- 
er  (1977)  which  is  based  on  Matsuno's  (1970.  1971)  formulation  for  ultra-long  waves. 

II.  MODEL 

A.  Equations 

The  approximate  equations  that  describe  the  symmetric  circulation  of  the  middle  atmos- 
phere can  be  derived  with  the  scaling  arguments  used  by  Leovy  (1964).  (cf.  Holton.  1975).  As 


our  equations  are  essentially  those  obtained  by  Leovy  (1964)  with  only  slight  differences  in  the 
final  formulation,  we  merely  summarize  the  salient  assumptions  and  arguments. 


The  equilibrium  atmosphere  is  assumed  to  be  motionless,  hydrostatic,  and  in  radiative 


balance  With  the  vertical  coordinate  ; * — 

P 


, where  is  a reference  pressure  taken  to  be 


lOOmb.  the  basic  state  equations  are 


a < <f>  > 
d: 


~ R < T> 


< H>  = < a > <T>  +<fi> 


(1) 


(2) 


where  < > denotes  a globally  averaged  quantity  (i.e.  zonally  and  meridionally  averaged),  a 
and  d are  the  Newtonian  cooling  coefficients,  Tis  the  temperature.  ^ is  the  geopoteniial,  R is 
the  dry  air  gas  constant,  and  ff  is  the  solar  heating. 

The  symmetric  circulation  is  described  by  meteorological  fields  denoted  by  an  overbar 
that  are  the  difference  between  the  zonally  averaged  fields  and  the  globally  averaged  fields. 
The  relevant  equations  are  (Leovy,  1964) 
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(4) 

(5) 

(6) 

(7) 


where  i7  and  vare  the  eastward  (westerly)  and  northward  (southerly)  wind  velocities.  0 is  lati- 
tude, a is  the  earth's  radius,  f)  is  the  earth's  rotation  rate,  k,  is  the  Rayleigh  friction  coefficient. 


3 


and  7 - — . the  'vertical  velocity'  The  last  term  on  the  RHS  of  Eqs.  <3)  and  (5)  represents 

the  divergences  of  the  zonally  averaged  eddy  momentum  and  heat  fluxes,  due  to  large  scale 
quasi-stationary  planetary  waves.  5 is  the  globally  averaged  static  stability  deflned  as 

R < T>  8 < r> 


S - 


or 


(8) 


H IS  the  diflerential  solar  heating  due  principally  to  the  eanh>sun  geometry.  The  term  T 
represents  diRerential  heating  or  cooling  due  to  O3  density  perturbations  that  accelerate  the 
thermal  relaxation  rate  in  the  upper  stratosphere  and  mesosphere  (Blake  and  Lindzen.  1773). 
An  effective  Newtonian  cooling  or  thermal  relaxation  rate.  is  defined  as  < a > a^. 

The  Eulenan  time  derivatives  in  Eqs.  (3)  • (S)  have  been  neglected  since  the  annual  oscillation 
frequency  associated  with  the  symmetric  circulation  is  much  less  than  the  damping  frequency. 
To  avoid  numerical  difficulties  which  arise  near  the  equator  the  term  k^v  hu  been  re- 
tained in  Eq.  (4). 

Eliminating  ii  and  v between  Eqs.  (3),  (4),  and  (7)  and  wfrom  (5)  and  (7),  we  arrive  at  a 
single  equation  for  <t 
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(9) 


where  y and  \ is  longitude.  The  eddy  flux  terms  are  discussed  in  Section  2E.  Eq.  (9) 

is  an  elliptic  differential  equation  and  may  be  solved  by  any  number  of  standard  numerical 
techniques.  We  used  a form  of  Crout  reduction  which  is  described  in  most  texts  on  numerical 
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methods  Adequate  convergence  was  obtained  with  vertical  grid  spacing  of  0 45  scale  heights 
and  horizontal  spacing  of  S latitude  degrees. 


B.  Boundary  Conditions 

The  boundary  condition  at  the  north  and  south  poles  is  obtained  Trom  Eq.  (4)  by  the  re- 
quirement that  u and  v must  vanish  there:  thus 


Along  the  lower  boundary  (;  - 0)  we  define  a geopotential  distribution  based  upon  the  obser- 
vational statistics  given  by  Taljaard  et  al.  (1968)  and  Labitzke  et  al.  (1972)  for  the  southern  and 
northern  hemispheres,  respectively.  Along  the  upper  boundary  (-~  125km)  we  assume  radia- 
tive equilibrium  is  valid  by  increasing  the  magnitude  of  to  include  thermal  conduction  as 
well  as  15m  CO;  radiation.  Thus  Eq.  (S)  reduces  to 


(11) 

at  the  upper  boundary.  Unfortunately  H cannot  be  precisely  determined  without  knowledge  of 
r since  the  local  densities  of  the  absorbing  constituents  O;  and  Oj  depend  on  T.  Accordingly 
we  assume  initially  that  7*  — 0,  compute  //,  and  solve  Eqs.  (9)  and  (11).  Using  the  new  T 
values,  we  update  H and  repeat  the  above  process  until  convergence  is  achieved.  In  practice 
only  two  or  three  iterations  are  required  to  obtain  0.2K  convergence  for  T. 

C.  Newtonian  Cooling 

The  Newtonian  cooling  coefficient,  — < a > 4-  a^.  is  a global  function,  while  the 
6 coefficient  is  defined  to  be  a function  of  height  only.  The  atmospheric  cooling  rate  due  to  IR 
emission  by  CO;  and  Oj  has  been  accurately  calculated  by  Dickinson  (1973a)  for  the 
30-100km  region.  Above  70km  the  accuracy  of  the  Newtonian  cooling  representation  breaks 


down  but  IS  still  retained  for  lack  of  a belter  representation.  We  estimate  the  asymptotic  value 


of  < II  > at  9Skm  to  be  the  COj  cooling  rate  in  its  fundamental  band  at  ISm  and  above  this 
level  by  the  thermal  conduction  relaxation  rate.  Below  JOkm  we  adopt  Trenbeth’s  (1973) 
values  for  < u > . The  resultant  < a > profile  is  displayed  in  Fig.  I.  The  ozone  photochemi* 
cal  acceleration  of  the  Newtonian  cooling  rate  is  based  on  Blake  and  Lindzen  (1973).  For  an 
atmosphere  with  odd  oxygen  chemistry  the  thermal  relaxation  rate  for  local,  small-scale  pertur- 
bations is  1400  ( <//>  + H)  (<T>  + f)  above  40km  where  the  recent  NBS  recom- 
mended rate  coefficients  have  been  adopted  (Hampson  and  Garvin.  1975).  In  the  real  atmos- 
phere where  odd  hydrogen  and  odd  nitrogen  chemistry  affects  ozone  concentrations  Blake  and 
Lindzen  concluded  that  the  above  value  of  should  be  reduced  by  approximately  one  half  to 


700 


( < r>  D' 


( < //>  4-  //) 


(12) 


Above  70km.  is  asymptotically  set  to  zero.  The  applicability  of  Eq.  (12)  to  global  scale  per- 
turbations. where  the  UV  transmissitivity  of  the  atmosphere  is  altered,  is  questionable  as  dis- 
cussed in  Section  3C. 


We  achieve  radiative  balance  in  the  atmosphere  as  represented  by  Eq.  (2)  in  the  follow- 
ing manner.  Since  the  globally  averaged  temperature  profile  is  fairly  well  known  and  the  m^or 
contributions  to  solar  heating,  <//>,  can  be  accurately  computed,  we  choose  <b>  in  Eq.  (2) 
to  obtain  the  required  radiative  balance  with  the  <a>  profile  shown  in  Fig.  1.  This  procedure 
is  more  accurate  than  any  attempt  now  to  calculate  <b>  from  first  principles  and  then  infer 
< T>  from  Eq.  (2).  We  have  found  that  the  derived  values  for  < b>  below  70km  are  in  ac- 
cord with  the  Dickinson  (1973a)  results. 

D.  Parameterization  of  the  Atmospheric  Heating  Rate 

In  order  to  numerically  simulate  the  middle  atmosphere  circulation,  accurate  and  compu- 
tationally efficient  parameterizations  of  solar  radiative  heating  by  O3  absorption  in  the  Hartley 
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region,  Huggins  bands,  and  (he  Chappius  bands  and  by  O;  absorption  in  (he  Schumann'Runge 
bands  and  continuum  are  required  A particularly  convenient  parameteri/aiion  was  introduced 
by  Lindzen  and  Will  11973)  who  noted  that  if  O)  absorption  cross  sections  are  either  constant 
or  have  an  exponential  variation  over  individual  wavelength  intervals,  the  atmospheric  heating 
rate  can  be  expressed  as  the  sum  of  simple  analytic  formulas  We  have  adopted  this  parame* 
terization  scheme  and  only  summarize  the  results  here,  details  may  be  found  in  Strobel  (1976). 

Below  approximately  80  km  the  0(^P)  formed  by  O]  and  Oj  dissociation  recombines 
quickly  and  the  actual  heating  rate  of  the  atmosphere  is  (he  total  photon  energy  absorbed 
Above  80  km  the  0(^P)  may  be  transported  significantly  before  recombining,  and  the  actual 
heating  rate  is  the  net  heating  rate  (defined  as  the  energy  absorbed  minus  the  chemical  energy 
stored  by  the  dissociated  products)  plus  the  chemical  energy  released  by  locally  recombining 
0(^P).  To  a good  first  approximation  we  can  estimate  the  actual  heating  rate  by  summing  the 
total  (net  + chemical)  Oj  and  the  net  heating  rates.  We  are  only  interest  tn  the  diflierential 
heating  rate  and  the  results  of  Park  and  London  (1974)  indicate  that  chemical  recombination 
heating  is  at  most  5K  day  above  80  km  and  that  transport  of  0(^P)  tends  to  minimize  the 
diflerential  chemical  recombination  heating  and  render  it  negligible  when  compared  to  the  net 
difierential  heating.  Our  heating  rates  agree  well  with  their  results  when  the  comparison  is 
made  at  the  same  density  levels. 

In  Table  1 the  total  heating  rates  for  Oj  and  the  net  heating  rates  for  Oj.  (?,.  are  tabulated 
by  spectral  region  in  units  of  ergs  cm  ~^s  and  are  related  to  Hby 

<H>  + H ~ — X Q, 

where  p is  the  mass  density  of  the  atmosphere.  The  column  densities  of  O2  and  Oj  (.Vj  and 
iVj.  respectively)  are  computed  by  an  8 point  Gaussian  quadrature  integration  scheme.  The 
number  densities  of  Oj  and  O3  are  denoted  by  rij  and  nj,  respectively. 
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The  time  scaJe  lor  motions  under  consideration  is  long  compared  to  a day  and  only  the 
diurnuUy  averaged  heating  rate  is  required.  We  have  adopted  approximation  9 4 of  Cogley  and 
Borucki  119761  to  accurately  and  quickly  evaluate  average  heating  rates.  We  have  obtained  an 
overall  accuracy  of  better  than  z 10  percent  when  compared  wun  'exact'  calculations. 

E.  Computation  of  Eddy  Fluxes 


The  eddy  heat  and  momentum  fluxes  in  equation  19)  are  computed  by  solving  a separate 
perturbation  equation  governing  planetary  wave  propagation,  similar  to  that  derived  by  Schoe- 
berl  and  Geller  (1977).  Stationary  planetary  waves  with  zonal  wavenumbers  one  and  two  are 
observed  to  have  sufficient  amplitude  in  the  stratosphere  at  mid  latitudes  to  account  for  greater 
than  90%  of  the  eddy  sutistics  (van  Loon  et  al.,  1973).  The  steady  state  wave  perturbations  are 
calculated  on  a grid  staggered  along  the  mid  points  of  the  mean  circulation  grid.  We  talte  the 
Rayleigh  friction,  (c,,  to  be  the  same  for  both  wave  and  mean  flow  geopoientials  but  Newtoni- 
an cooling  for  the  wave  is  approximated  by  <a>  since  planetary  waves  are  important  only  in 
the  winter  hemisphere  where  is  small.  The  planetary  wave  equation  is: 


(ifufn  + k,) 


•h  (itu/fl  + < a > ) 
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The  boundary  conditions  on  equation  (13)  are: 

It'  -0.  r - 16  (-125  km) 

4)'  “/■(*♦).  z 0 (—16km) 

where  Fin)  is  the  climatological  forcing  given  by  van  Loon  et  al.  (1973). 

The  steady  state  contribution  of  planetary  wave  eddy  fluxes  on  the  mean  circulation  is 
evaluated  iteratively  by  computation  of  the  wave  amplitude  in  the  presence  of  the  zonal  flow  at 
that  iteration  followed  by  the  adjustment  of  the  mean  circulation  to  the  divergence  of  the  eddy 
fluxes  just  computed.  This  iteration  procedure  is  continued  until  both  the  wave  amplitude  and 
the  mean  circulation  adjustments  decreas  .o  less  than  5%. 

F.  Model  Atmosphere 

CIRA  (1972)  was  the  model  atmosphere  adopted  for  < T>  and  the  O and  O2  mixing  ra- 
tios. The  O3  density  distribution  is  based  on  Krueger  (1973).  Other  Oj  density  distributions 
were  considered  above  75  km  and  will  be  discussed  below  in  Section  3F  The  latitudinal  varia- 
tion of  O3  density  along  constant  pressure  surfaces  has  been  neglected  in  the  present  model 
since  differential  solar  heating  depends  primarily  upon  earth-sun  geometry. 

III.  THE  SOLSTITIAL  CIRCULATION 

A principal  goal  of  this  study  is  to  assess  the  quantitative  effects  of  various  physical 
processes  on  the  mean  circulation  of  the  middle  atmosphere.  We  proceed  by  comparing  calcu- 
lations where  these  processes  have  been  added  or  deleted  one  at  a time  to  see  how  the  total 
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circulaiion  is  altered.  The  most  sensitive  indicator  of  the  mean  circulation  is  the  zonal  wind  at 


solstice.  In  each  test  calculation  described  in  this  section  we  show  only  the  mean  zonal  wind 
profile  but  give  other  diagnostic  variables  at  selected  altitudes  m Table  11.  All  diagnostic  vari- 
ables are  shown  for  the  final  calculations. 

A.  The  Heating  Rale 

The  perturbation  heating  rate.  H,  at  solstice  is  shown  in  Figure  2,  while  the  globally  aver- 
aged rate.  <H>.  is  illustrated  in  Figure  3.  In  the  stratosphere  O3  absorption  dominates  the 
heating  rate  whereas  absorption  is  responsible  for  the  large  heating  rate  above  90  km.  Note 
that  to  first  order.  H is  antisymmetric  about  the  equator.  There  is.  however,  a non-negligible 
symmetric  component  of  H. 

^ Leovy  (1964)  showed  that  the  antisymmetric  component  of  the  heating  drives  a single, 

direct  pole  to  pole  circulation  cell  with  upward  motion  in  the  summer  hemisphere  and  down- 
ward motion  in  the  winter  hemisphere.  Conversely,  symmetric  heating  forces  two  cells  with 
upward  motion  over  equatorial  regions  and  downward  motion  at  the  pole.  Coriolis  torques  act 
on  the  meridional  circulation  to  produce  strong  zonal  winds.  For  antisymmetric  heating  the 
zonal  wind  jets  are  easterly  in  the  summer  hemisphere  and  westerly  in  the  winter  hemisphere, 
whereas  symmetric  heating  produces  westerly  jets  in  both  hemispheres.  The  heating  shown  in 
Figure  2 will  yield  features  of  both  components. 

B.  Rayleigh  Friction 

The  principal  adjustable  coefficient  in  this  model  is  the  Rayleigh  fnciion  coefficient,  k,. 
There  is.  unfortunately,  no  self  consistent  theory  which  can  be  used  to  estimate  its  value.  Leo- 
vy (1964)  adopted  a value  of  k,~  0.1  day  which  Green  (1972)  criticized  as  being  too  large. 
To  illustrate  the  sensitivity  of  the  zonal  Jets  to  the  Rayleigh  friction  we  present  calculations  for 
three  constant  values  of  k,"' ; 30,  15  and  7.5  days  in  Figure  4a,  b.  c,  respectively. 
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In  the  absence  ot'  RayleiKh  friction,  the  meridional  circulation  vanishes  if  the  vertical 
velocity  IS  zero  at  the  lower  boundary,  as  is  the  case  here,  and  the  temperature  held  is  strictly 
in  radiative  equilibrium.  Under  these  conditions  a large  meridional  temperature  gradient  is 
produced  which  generates  a vertical  shear  in  the  zonal  wind,  consistent  with  the  thermal  wind 
equation.  This  case  is  illustrated  m Fig.  4a  where  the  strength  of  the  jets  increases  with  alii' 
lude  up  to  ~ 1 10  km.  Since  is  small,  but  finite  the  wind  contours  close  at  high  altitudes  in- 
dicating the  presence  of  a weak  pole  to  pole  circulation.  The  zonal  wind  structure  shown  in 
Fig.  4a  bears  little  resemblance  to  observations  reported  in  .Vfurgatroyd  (1969)  and  CIRA 
(1965.  1972). 

For  larger  values  of  k,  a stronger  vertical  field  develops  in  response  to  the  external  heat- 
ing in  the  stratosphere  and  reduces  the  meridional  temperature  graidient  by  adiabatic  processes. 
More  specifically,  upward  motion  at  the  summer  pole  decreases  the  temperature  from  its  radia- 
tive equilibrium  value  and  downward  motion  over  the  winter  pole  raises  the  temperature  above 
its  radiative  equilibrium  value.  Consequently  the  stratospheric  zonal  wind  shear  weakens.  In 
the  mesosphere  the  strength  of  the  vertical  motions  are  sufficiently  large  to  reverse  the  meri- 
dional temperature  gradient  imposed  by  radiative  balance  and  thus  close  the  zonal  wind  con- 
tours below  the  mesopause.  In  Figs.  4b  and  4c  this  result  is  demonstrated.  The  jet  maxima  oc- 
cur between  60  and  80  km  in  agreement  with  observations,  although  the  jet  intensities  are  too 
great  in  Fig.  4b  with  ■ 15  days. 

The  temperature  gradient  at  the  mesopause  as  indicated  in  Table  2 is  a good  indicator  of 
the  magnitude  of  the  vertical  motion  field.  Another  important  feature  which  distinguishes 
these  three  calculations  is  the  separation  of  the  jet  maxima  in  each  hemisphere  into  two  lobes 
with  one  peak  in  the  mesosphere  and  the  other  in  the  lower  thermosphere.  Associated  with 
these  lobes  are  the  two  major  heating  regimes;  O3  and  insolation,  respectively  CIRA 
(1965,  1972)  suggest  easterly  winds  above  the  mesopause  in  equatorial  regions  (where  our 
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quasigeostrophic  model  is  inadequate)  and  westerlies  at  mid  latitudes  in  both  hemispheres  Ob- 
servations for  this  region  are  difficult  to  obtain  because  of  atmospheric  tidal  motions.  Calcula- 
tions by  Dickinson  et  al.  (1977)  of  the  thermospheric  circulation  also  indicate  the  presence  of 
jets  similar  to  those  shown  in  Figure  4c  with  amplitudes  of  80-120  ms  peaking  at  120  km 
and  about  60°  latitude. 

From  the  results  shown  in  Figure  4 it  is  clear  that  large  values  of  > 0.1  day  are  re- 
quired above  the  stratopause  to  produce  the  observed  wind  and  temperature  structure.  For  the 
lower  stratosphere.  Lindzen  (private  communication,  1977)  has  noted  that  k~^  values  of  10  to 
IS  days  would  severely  dampen  the  quasi-biennial  oscillation  (Holton  and  Lindzen,  1971). 
These  two  constraints  suggest  a Rayleigh  friction  profile  with  the  following  functional  form: 

(14) 

where  c — 4.43  and  kg”'  “ 40  days.  The  computed  zonal  wind  structure  with  this  Rayleigh 
friction  profile  is  shown  in  Figure  S.  The  agreement  between  the  observed  zonal  wind  struc- 
ture and  this  case  is  better  than  was  achieved  with  a constant  k,  value  (cf.  Figure  4c).  Note 
also  from  Table  2 that  the  computed  mesopause  temperatures  are  also  in  better  agreement  with 
the  observations.  The  warmer  temperatures  in  the  winter  mesosphere  and  cooler  temperatures 
in  the  summer  mesosphere  are  associated  with  a more  vigorous  meridional  circulation  above 
the  stratopause  as  a result  of  the  enhanced  Rayleigh  friction.  Vertical  motion  fields  generated 
by  heating  at  the  stratopause  increase  in  intensity  with  altitude  because  k,  and  ware  small  and 
cannot  "short  circuit"  the  w field  by  continuity  considerations  in  the  stratosphere.  Other  func- 
tional forms  for  Rayleigh  friction  are  not  ruled  out  by  Equation  (14).  We  find  only  two  criteria 
to  be  important  in  generating  realistic  mean  fields  in  the  middle  atmosphere:  (a)  that 
k,“'  > 40  days  in  the  lower  stratosphere  and  (b)  k,~'  ~-  10  days  above  the  stratopause. 


12 


C.  Photochemical  Damping 


We  now  include  photochemical  acceleration  of  Newtonian  cooling  with  the  rate  given  by 
equation  (12)  and  retain  the  exponential  k,  profile.  The  results  of  this  calculation  are  shown  in 
Figure  6.  From  a comparison  of  Figures  5 and  6 it  is  apparent  that  the  easterly  jet  in  the  sum- 
mer stratosphere  is  more  strongly  affected  with  a 15  ms  reduction  at  the  peak.  Additionally, 
the  summer  stratopause  temperature  is  lowered  by  about  lOK.  Photochemical  acceleration 
effectively  reduces  the  meridional  temperature  gradient  in  the  summer  hemisphere  since  is 
proportional  to  < H>  + H.  As  expected,  the  zonal  jet  in  the  summer  mesosphere  is  now  less 
intense  than  observed  which  may  indicate  that  smaller  values  of  k,  are  appropriate  for  that  re- 
gion or  that  has  been  overestimated  by  equation  (12). 

The  rale  for  was  derived  for  local,  small  scale  perturbations  where  the  UV  transmis- 
sivity of  the  atmosphere  is  independent  of  the  perturbation  (cf.  Blake  and  Lindzen,  1973).  For 
global  scale  perturbations  with  vertical  scales  greater  than  the  atmospheric  scale  height,  the  UV 
transmissivity  is  significantly  affected  for  optical  depths  > 1.  Estimates  for  the  latter  case  indi- 
cate that  at  unit  optical  depth,  the  perturbation  to  the  heating  rate  is  zero  (i.e.  ^-0)  because 

(he  change  in  the  heating  rate  due  to  local  Oj  density  penurbations  is  precisely  cancelled  by 
the  change  in  UV  transmissivity  (Strobe),  1977).  Thus  for  global  scale  perturbations  the  rate 
given  by  Eq.  (12)  is  appropriate  only  for  the  optically  thin  case  where  the  heating  rate  and 
hence  are  small.  As  a consequence  we  have  neglected  the  photochemical  acceleration  of 
Newtonian  cooling  in  subsequent  calculations. 

D.  Realistic  Lower  Boundary  Conditions 

In  our  previous  calculations  the  lower  boundary  condition  was  6 in.  0)  ->0.  We  now  in- 
clude a geopotential  perturbtion  defined  by  observational  statistics  compiled  by  Taljaard  et  al. 
(1969)  and  Labitzke  et  al.  (1972).  Figure  7 presents  calculations  for  m,  v,  w,  and  T when  sum- 
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mer  is  in  the  Southern  hemisphere,  while  Figure  8 is  for  winter.  The  u and  w helds  exhibit 
the  most  noticable  asymmetries  and  thus  the  5 and  T fields  are  not  presented  for  the  later  cal- 
culation We  have  used  the  exponential  Rayleigh  friction  profile  but  have  neglected  photo- 
chemical acceleration  of  Newtonian  cooling. 

It  is  apparent  from  these  figures  that  the  tropopause  height  field  induces  considerable 
hemispherical  asymmetry  into  the  stratosphere  and  mesosphere.  Particularly  remarkable  is  the 
upward  extension  of  tropospheric  westerlies  into  the  stratospheric  easterly  jet  in  Figure  7a.  No 
such  extension  occurs  in  the  equivalent  calculations  for  summer  in  the  Northern  hemisphere 
(cf.  Figure  8a).  Overall,  realistic  lower  boundary  conditions  have  a major  effect  only  below  50 
km  where  t^e  westerlies  are  intensified  in  both  winter  hemispheres,  the  summer  easterlies  in 
the  Southern  hemisphere  are  weakened,  and  the  Northern  hemisphere  summer  easterlies  are 
intensified. 

The  relatively  weak  influence  at  the  boundary  forcing  on  the  mean  circulation  is  con- 
sistent with  the  analytic  solutions  of  equation  (9)  which  have  been  studied  by  Leovy  (1964) 
and  Longuet-Higgins  (1968).  Longuet-Higgins  denotes  these  solutions  as  Type  7 solutions  to 
Laplace's  tidal  equation.  The  equivalent  depths  associated  with  the  solutions  are  small  and 
negative  so  boundary  effects  will  decay  rapidly  with  height  as  our  results  indicate. 

E.  Eddy  Fluxes 

We  now  consider  the  influence  of  vertically  propagating  stationary  planetary  waves  with 
zonal  wave  numbers  m ■■  1 and  /n  - 2.  The  vertical  structure  of  these  disturbances  is  con- 
trolled by  the  meridional  gradient  of  the  potential  vorticity,  and  the  strength  and  the 

fjtr 

shape  of  the  polar  night  jet  (Matsuno,  1970;  Schoeberl  and  Geller.  1977),  The  function 
cos  » is  shown  in  Figure  9 for  the  zonal  wind  structure  presented  in  Figure  7.  where  winter 
IS  m the  Northern  hemisphere 
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A chungc  in  the  sign  of  is  a necessary  condition  Cor  wave  instability  (Charney  and 

Stern,  l')62)  The  negative  region  above  80  km  in  the  winter  mesosphere  on  the  right  hand 
side  of  Figure  results  from  negative  shear  above  the  jet  maximum  as  discussed  previously  by 
Dickinson  (1973b)  and  Simmons  (1974)  They  suggest  planetary  waves  may  be  baroclinically 
unstable  m this  region.  Simmons  (1977)  has  also  examined  the  possibility  of  such  instability 
in  the  summer  mesosphere  as  a result  of  positive  shear  in  the  stratospheric  easterlies.  The 
growth  of  unstable  waves  in  either  region  is  opposed  by  the  large  radiative  and  photochemical 
damping  rales  so  no  Rrm  conclusion  has  been  reached  on  the  importance  of  unstable  waves  in 
the  mesosphere.  In  order  to  compute  the  influence  of  vertically  propagating  planetary  waves  in 
steady  state  equilibrium  with  the  mean  circulation  we  suppress  instabilities  by  giving 
cos  H a positive  value  of  10  s where  a sign  change  occurs. 

For  an  arbitrarily  large  wave  amplitude  forced  at  100  mb  there  is  no  guarantee  that  the 
Iterative  scheme  will  converge.  The  modification  of  the  mean  flow  by  wave  flux  convergences 
may  enhance  these  fluxes  at  the  next  iteration  and  thus  lead  to  a numerically  non -convergent 
sequence.  For  example,  if  easterlies  appear  in  the  region  dominated  by  the  wave  field,  conver- 
gence will  be  unlikely  due  to  critical  level  effects  (Dickinson.  1969;  Matsuno.  1971).  Note  also 
that  we  have  not  included  the  displacement  of  the  equatorial  easterlies  by  the  quasi  biennial 
oscillation  so  that  possible  critical  level  effects  have  been  neglected  at  low  latitudes. 

For  m — 1 we  could  not  obtain  numerical  convergence  using  climatological  wave  ampli- 
tudes given  by  van  Loon  et  al.  (1973).  After  the  first  iteration,  easterly  winds  appeared  near 
the  winter  pole  mesopause  and  the  iterative  procedure  diverged  due  to  the  build  up  of  eddy 
heal  and  momentum  flux  convergence  at  the  critical  level.  Numerical  convergence  was  ob- 
tained when  the  wave  amplitude  was  uniformly  reduced  to  1/3  of  the  climatological  value  at 
the  lower  boundary  to  a peak  amplitude  of  45  gpm.  With  reduced  wave  amplitude,  the  »i  “ 1 
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eddy  excried  a negligible  mlluencc  on  the  mean  circulaiion  since  iis  flux  convergence  is  pro- 
portional to  i<b  To  illustrate  the  impact  of  the  /»/  — 1 eddy  with  the  full  climatological  am- 
plitude at  100  mb  we  show  the  meteorological  fields  after  the  first  iteration  when  the  eddy 
fluxes  were  added  as  heat  and  momentum  sources  in  Fig.  10.  We  have  found  that  the  conver- 
gence of  northward  eddy  heat  fluxes  is  strongest  at  about  85  km  and  poleward  of  TS'N  with 
secondary  convergence  occurring  at  30  km  poleward  of  6S°N  associated  with  local  maxima  of 
wave  amplitude  (Figure  lOe).  The  convergence  of  zonal  momentum  flux  occurs  at  approxi- 
mately the  same  altitudes  and  latitudes  with  an  additional  divergence  occuring  at  65  km  pole- 
ward  of  65°N.  The  primary  eddy  forcing  of  the  mean  flow  appears  to  be  the  northward  heat 
flux  convergence  which  sets  up  a local  circulation  cell  near  the  winter  pole  (Figure  10c,  d)  and 
generates  the  strong  easterlies  there  (Figure  10a).  The  intensity  of  the  polar  night  jet  is  re- 
duced by  ~10%  as  a result  of  eddy  heat  and  momentum  flux  convergences,  with  the  heat  flux 
convergence  at  30  km  being  the  dominant  contributor,  (n  the  winter  polar  mesosphere  the 
temperature  has  increased  by  — 70K  due  to  the  m — 1 eddy.  This  temperature  extreme  is,  of 
course,  not  observed  and  the  circulation  pattern  shown  in  the  upper  parts  of  Figure  10  should 
be  regarded  as  illustrative.  The  points  demonstrated  by  this  calculations  are  (a)  a reduction  of 
the  jet  strength  by  stratospheric  eddy  heat  convergence  and  (b)  the  appearance  of  easterly 
winds  in  the  upper  mesosphere  due  to  eddy  heat  convergence  in  the  mesosphere.  This  later 
region  might  be  regarded  as  an  initiation  point  for  critical  levels  which  could  later  progress 
downward  to  generate  sudden  stratospheric  warmings  (Matsuno,  1971;  Holton,  1976). 

The  penetration  of  the  polar  night  jet  by  the  m — 1 planetary  wave  shown  in  Figure  lOe 
is  consistent  with  calculations  made  by  Schoeberl  and  Geller  (1977)  who  showed  that  the  m - 
1 wave  IS  generally  not  trapped  except  for  very  strong  jets.  The  numerical  divergence  of  the 
Iteration  method  may  also  indicate  that  the  stratospheric  m — 1 wave  never  attains  a steady 
solution  but  will  vacillate  in  amplitude  throughout  the  winter.  Vacillation  behavior  of  plane- 
tary waves  has  been  examined  by  Holton  and  Mass  (1976).  The  reduced  amplitude  which  we 
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found  would  give  convergence  may  be  related  to  the  minimum  amplitude  required  for  such 
stratospheric  vacillations,  however,  we  note  that  our  minimum  amplitude  is  about  three  times 
smaller  than  that  deduced  by  Holton  and  Mass.  This  discrepancy  is  probably  due  to  the  fact 
that  the  winter  mean  zonal  wind  profile  used  by  Holton  and  .Mass  does  not  decrease  above  60 
km. 

For  planetary  wave  m — 2 convergence  could  be  obtained  for  climatological  forcing.  The 
wave  amplitude  for  m -•  2 is  shown  in  Figure  11;  note  that  this  wave  is  trapped  in  the  lower 
stratosphere  with  a smaller  amplitude  than  m « I (Figure  lOe).  The  effect  of  m ••  2 on  the 
mean  circulation  is  negligible  as  might  be  expected  from  the  weaker  wave  amplitude.  The 
trapping  of  this  wave  by  the  polar  night  jet  is  also  consistent  with  calculations  made  by 
.Vlatsuno  (1970)  and  Schoeberl  and  Geller  (1977).  Clearly,  for  m “ 2.  the  eddy  heat  flux  con- 
vergence will  be  largest  in  the  lower  stratosphere,  but  the  mean  circulation  associated  with  this 
convergence  is  so  weak  that  less  than  a 1%  reduction  in  the  jet  amplitude  occurs. 

F.  The  Local  Static  Stability  as  a Function  of  Ozone  Density  ProBles 

In  the  course  of  this  investigation  we  discovered  that  the  radiative  heating  associated  with 
certain  ozone  profiles  generated  perturbation  equilibrium  temperatures  sufficiently  large  in  the 
polar  mesopause  region  to  drive  the  local  static  stability  negative  This  is  a necessary  condition 
for  convection  and  invalidates  the  use  of  a globally  averaged  static  stability  in  cur  model.  Fig- 
ure 12  shows  three  ozone  profiles  tested  by  our  model:  Profile  1 based  on  Hays  and  Roble 
(1973);  Profile  II  extrapolated  from  Krueger  (1973)  and  used  in  our  previous  calculations;  and 
Profile  III.  an  intermediate  profile.  Profile  1 yields  locally  large  heating  regions  —2  scale 
heights  below  the  Schumann-Runge  band  heating  region.  The  rapid  variation  of  <(«>  and  its 
small  value  combined  with  the  large  altitude  variation  of  H produce  a radiative  equilibrium 
temperature  profile  with  a super-adiabatic  lapse  rate.  Such  variations  are  also  present  for 
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Proliles  II  and  III  but  wiih  smaller  amplitude  and  generate  only  sub-adiabatic  lapse  rates.  We 
note  that  if  Profile  I were  produced  by  ozone  photochemistry  and  our  choice  of  the  Newtonian 
cooling  rate  is  appropriate,  then  the  convection  and  associated  turbulence  generated  by  this 
large-scale  thermal  forcing  may  be  self-sustaining  since  the  chemical  and  convection  time  con- 
stants are  comparable  ( — 10  min).  As  a result  of  this  convection,  the  summer  polar  mesopause 
temperature  could  be  lowered  sufficiently  to  condense  water  vapor  and  form  noctilucent  clouds. 
The  convection  would  also  generate  the  vigorous  vertical  motion  required  by  Reid  11975)  in 
his  model  of  these  ice  clouds. 

IV.  THE  MEAN  CIRCULATION  AT  EQUINOX 

Using  the  Rayleigh  friction  profile  given  by  Eq.  (14).  lower  boundary  geopotential  fields 
from  Taljaard  et  al.  (1969)  and  Labitzke  et  al.  (1972)  for  April  conditions,  and  *0.  we 
have  computed  the  mean  circulation  at  equinox  which  is  shown  in  Fig.  13.  The  perturbation 
heating  rate  (Fig.  13a)  is  symmetric  about  the  equator  as  expected  with  peak  heating  rate  of 
—3  K day  at  the  equator  and  maximum  cooling  rate  of  — IIK  day  at  the  poles  in  the 
ozone  layer.  The  slight  asymmetry  in  the  zonal  jets  results  from  the  asymmetrical  lower  boun- 
dary conditions.  Both  the  temperature  and  zonal  wind  profiles  are  in  good  agreement  with  ob- 
servations (CIRA.  1965,  1972);  however.  Table  2 indicates  that  the  calculated  jets  occur  at  a 
lower  altitude.  The  meridional  circulation  pattern  illustrates  the  double  cell  structure  charac- 
teristic of  the  symmetric  heating  at  equinox:  rising  motion  over  the  equatorial  region  and  sink- 
ing motion  at  the  poles.  We  conclude  that  the  extratropical  mean  circulation  at  equinox  and 
solstice  are  driven  by  the  same  physical  processes.  Equatorial  motions  are  not  adequately 
described  by  our  quasigeostrophic  model  and  no  attempt  is  made  to  compare  computed  equa- 
torial motions  with  observations 

At  midlatitudes.  Dickinson  (1969)  has  shown  that  vertically  propagating  planetary  waves 
undergo  strong  radiative  damping  at  equinox  and  the  extratropical  effects  of  planetary  wave 
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eddy  stresses  are  regarded  as  negligible.  Rayleigh  friction  may  thus  be  weaker  during  eguinos 
which  would  raise  the  location  of  the  calculated  ?onal  jets  and  remove  the  discrepancy  noted 
previously 

V.  SUMMARY 

Using  a numerical  quasigeostrophic  model  we  conclude  that  the  essential  features  of  the 
2onally  averaged,  extratropical  circulation  can  be  adequately  understood  if  Dickinson's  (1973a) 
Newtonian  cooling  rate  is  the  appropriate  thermal  relaxation  rate  and  friction  is  described  by 
the  Rayleigh  coefficient  given  in  Eq.  (14).  In  particular  we  find  the  direct,  radiatively  driven 
pole  to  pole  circulation  at  solstice  is  sufficient  to  account  for  the  cold  summer  mesopause  and 
the  warm  isothermal  winter  mesosphere  with  associated  zonal  jets  of  realistic  magnitude,  thus 
confirming  Leovy’s  (1964)  original  hypothesis. 

The  climatological  heat  and  momentum  fluxes  associated  with  planetary  wave  nugiber  2 
have  a negligible  effect  on  the  mean  circulation.  For  planetary  wave  in  - 1 we  could  not  ob- 
tain a steady  stale  contribution  of  eddy  forcing  to  the  mean  circulation  given  the  observed 
climatological  wave  amplitude  at  the  tropopause.  Instead  we  suggest  that  the  stratospheric 
response  to  this  magnitude  of  m “ 1 forcing  will  be  a .vacillation  behavior  as  previously  dis- 
cussed by  Holton  and  Mass  (1976).  According  to  our  calculations  the  threshold  amplitude  for 
this  response  is  only  1/3  of  the  climatological  forcing,  a factor  of  3 less  than  the  threshold  ob- 
tained by  Holton  and  Mass.  Below  the  threshold  value  the  m - 1 planetary  wave  has  a negligi- 
ble effect  on  the  mean  circulation  of  the  middle  atmosphere. 

Our  computations  suggest  that  the  thermal  damping  rate  in  the  vicinity  of  the  stratopause 
cannot  substantially  exceed  Dickinson's  Newtonian  cooling  rate.  The  photochemical  accelera- 
tion of  the  thermal  damping  rate  calculated  by  Blake  and  Lindzen  (1973)  is  inapplicable  to  the 
global,  large  vertical  scale  motion  considered  here  because  LiV  transmissivity  effects  are  impor- 
tant. 
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Another  interesting  result  of  our  calculations  is  the  inference  that  secondary  peaks  in  the 
O3  density  near  the  mcsopause  will  thermally  force  the  summer  polar  mesopause  region  to  be 
convectiveiy  unstable  and  signilicantly  lower  its  temperature.  This  increases  the  possibility  of 
water  vapor  condensation  to  form  ice  clouds.  In  addition  convection  will  supply  the  vigorous 
vertical  motion  required  by  Reid  (1975)  in  his  model  of  noctilucent  clouds. 

ACKNOWLEDGEMENTS 

The  authors  would  like  to  thank  Professor  R S.  Lindzen  for  useful  comments  and  Harry 
S.  van  Loon  for  providing  planetary  wave  data.  The  major  part  of  this  research  was  supported 
by  the  Naval  Air  Systems  Command  with  additional  support  from  the  Upper  Atmospheric 
Research  Office  of  the  National  Aeronautics  and  Space  Administration. 

REFERENCES 

Blake,  D.  and  RS.  Lindzen,  1973:  Effect  of  photochemical  models  on  calculated  equilibria  and 
cooling  rates  in  the  stratosphere,  Mon.  IVea.  Rev.,  101,  783-802. 

Charney,  J.G.  and  M.E.  Stern,  1962:  On  the  stability  of  internal  baroclinic  jets  in  a rotating  at- 
mosphere. J.  Atmos.  Set.,  19,  159-172. 

CIRA  1965:  COSPAR  International  Reference  Atmosphere,  Amsterdam,  North-Holland,  313  pp. 
CIRA  1972:  COSPAR  International  Reference  Atmosphere,  Berlin,  Akedemie-Verlag,  450  pp. 
Cogley,  A.C.  and  WJ.  Borucki,  1976:  Exponential  Approximation  for  Daily  Average  Solar  Heat- 
ing or  Photolysis.  J.  Atmos.  Sci.,  33,  1347-1356. 

Dickinson,  R.E.,  1969:  Vertical  propagation  of  planetary  Rossby  waves  through  an  atmosphere 
with  Newtonian  cooling.  J.  Geophys.  Res.,  74,  929-938. 

Dickinson,  R.E..  1973a:  Method  of  parameterization  for  infrared  cooling  between  altitudes  of 
30  and  70  kilometers,  J.  Geophys.  Res.,  78.  4451-4457. 

Dickinson.  R.E.,  1973b:  Baroclinic  instability  of  an  unbounded  zonal  shear  flow  in  a compressi- 
ble atmosphere.  J.  Atmos.  Sci.,  30,  1520-1527. 


20 


Dickinson.  R E..  E C Ridley,  and  R.G,  Roble.  1977:  Meridional  Circulation  in  ihe  Thermos- 
phere II  Solstice  conditions  J.  Aumn.  SV/..  34.  178-192. 

Green.  J.S.A..  1972:  Large  Scale  motion  in  the  upper  stratosphere  and  mesosphere;  an  evalua- 
tion of  data  and  theories.  Phil.  Tram.  Roy.  Sik.  Lout/..  271.  577-583. 

Hampson.  R.F..  Jr.,  and  D Garvin.  1975;  Chemical  kinetic  and  photochemical  data  for  model- 
ling atmospheric  chemistry  NBS  Tech.  Note.  866.  1 13  pp. 

Hays.  P.B.  and  R.G  Roble,  1973;  Observation  of  mesospheric  ozone  at  low  latitudes.  Planet. 
Space  Sci..  21,  273-279. 

Holton,  J.R.,  1975;  The  Dynamic  .VteieoroloKy  of  the  Stratosphere  and  .Sfesosphere  .Vfeteor.  .Won..  15. 
Boston,  Amer.  Met.  Soc..  216  pp. 

Holton.  J.R..  1976;  A semi-spectral  numerical  model  for  wave-mean  flow  interactions  in  the 
stratosphere;  applications  to  sudden  stratospheric  warnings.  J.  Atmos.  Sci.,  33.  1639-1649. 
Holton,  J.R.  and  R5.  Lindzen,  1972:  An  updated  theory  for  the  quasi-biennial  cycle  of  the 

C 

tropical  stratosphere.  J.  Atmos.  Sci.,  32,  712-719. 

Holton.  J.R.  and  C.  Mass.  1976.  Stratospheric  vacillation  cycles.  J.  Atmos.  Set.,  33,  2218-2225. 
Krueger,  AJ..  1973:  The  mean  ozone  distribution  from  several  series  of  rocket  soundings  to  52 
km  at  latitudes  from  58°S  to  64'’N.  Pure  and  Appl.  Geophys.,  106-108.  1272-1280. 

Labitzke.  K.  and  Collaborators.  1972;  Climatology  of  the  Stratosphere  in  the  northern  hemisphere 
Part  I,  Heights,  temperatures  and  geostrophic  resultant  wind  speeds  at  100.  50,  30.  and  10  mb. 
.Weteorol.  Abhaitdl.  100,  6 pp. 

Leovy,  C.B.,  1964;  Simple  models  of  thermally  driven  mesospheric  circulations.  J.  Atmos.  Sci, 
21,  327-341. 

Lindzen,  R.S.  and  D.I.  Will.  1973;  An  analytic  formula  for  heating  due  to  ozone  absorption,  J. 
Atmos.  Sci..  30,  513-515. 

1 


Longuc(-Higgins.  M S.,  1968:  The  eigenluncoons  of  Laplaces'  tidal  equations  over  a sphere 
Phil.  Tram.  Roy.  Soi.  LonJoii.  A262.  51 1-607, 

Matsuno,  T..  1971:  A dynamic  model  of  the  stratospheric  sudden  warming.  J.  .4imos.  Sii..  28, 
1479-1492. 

Matsuno.  T . 1970  The  vertical  propagation  of  stationary  planetary  waves  in  the  winter  north- 
ern hemisphere.  / .4iinos.  Sci.,  27,  871-883. 

Murgatroyd,  RJ..  1969:  The  structure  and  dynamics  of  the  stratosphere  The  Global  emulation 
ol  the  .Atmosphere.  G.A.  Corby,  Ed.,  London,  Roy.  Meteor.  Soc.,  159-195. 

Park,  J.H.  and  J.  London,  1974:  Ozone  photo  chemistry  and  radiative  heating  of  the  middle  at- 
mosphere. J.  Atmos.  Svi..  31,  1898-1961. 

Reid.G.C.,  1975:  Ice  clouds  at  the  summer  polar  mesopause.  J.  Atmos.  Sci..  32.  523-535. 

Schoeberl,  M.R.  and  M.A.  Geller.  1977:  A calculation  of  the  structure  of  stationary  planetary 
waves  in  winter,  J.  Atmos.  Sci..  34,  in  press. 

Simmons.  AJ..  1974;  Baroclinic  instability  at  the  winter  stratopause  Quart.  J.  Roy.  Meteor.  Soc. 
100,  531-540. 

Simmons.  A.J.,  1977;  Baroclinic  instability  in  the  summer  mesosphere.  Quart.  J.  Roy.  Meteor. 
Soc..  103.  211-215. 

Strobel,  D.F..  1976;  Parameterization  of  the  atmospheric  heating  rate  from  15  to  120  km  due  to 
Oj  and  O3  absorption  of  solar  radiation.  NRL  Memo  Rept.  3398,  Naval  Research  Labora- 
tory. Washington  DC. 

Strobel,  D.F..  1977:  Photochemical-Radiative  Damping  and  Instability  in  the  Stratosphere,  Geo- 
phys.  Res.  Lett.,  submitted. 

Taljaard.  J.J.,  H.  van  Loon,  H.L.  Crutcher  and  R.L.  Jenne,  Climate  oT  the  upper  air.  Part  I • 
Southern  hemisphere,  NAVAIR  50-IC-55,  3 pp. 


Trenbeih.  K E . l'J7J  Global  model  of  the  general  circulation  of  the  atmosphere  below  75  km 
with  animal  heating  cycle  Miw  Ufa  Rev..  101,  306-322. 
van  Loon,  H.,  R.L.  Jenne,  and  K Labitzke  1973  Zonal  harmonic  standing  waves,  J.  Geophvs. 
Rf',..  78.  4463-4471 


23 


Table  I - Parameterization  of  the  Atmospheric  Heatinf  Rate 


Spectral  Region 
Chappius  bands 

Huggins  bands 

Hartley  region 
Herzberg  continuum 

Schumann-Runge  bands 

Schumann-Runge  continuum 

*0,  is  in  ergs  cm  s n,  in 


Parameterization* 

— -1.05  X 10 exp( -2.85  x 10 
"3 

— -^14.66  X 10  ^-7.8  X 10  ^xp(-1.77  x 10-'’Vj) 

"3  "3 

-3.88  X 10^  exp  ( -4.22  x 10  iVj  )I 

— - 4.8  X 10  -‘♦exp  { -8.8  x 10  N-,  ) 

"3 

Qh.  - (9.9  x 10  /I2  + 7 J5  X 10  flj  ) X 

!xp  ( -6.6  X 10  /Vj  -4.9  X 10  A'j  ) 

1 


" 0.67  ACj  + 3.44  X 10’ 
< 10'*  c/n^, 

- - 2.43  X 10 


exp 


- •-zr  10-98  exp  ( -2.9  x 10  N-, ) -0^5 

Nj 

1.7  X 10-'*Nj)  -0.43  exp  ( -1.15  x 10  Nj  )) 
‘.'A,  in  cm 
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IKm) 


OTITUOe  lOfGl 


Fi|.  4 - Th«  tiMan  zonal  wind  dUtribution  computed  uiinf  constant  Raykiyh  fne- 
tion  coalTicimls.  Fof  • <al  30  days,  tb)  15  days.  (c(  days.  Wind  velocity 
It  in  ms**. 


28 


LATITUDE  (OEGI 


5 - The  mean  zonal  wind  distribution  computed  usintt 
the  Rayleigh  friction  profile  given  by  equation  (14). 


•7  5 58  29  0 29  58  87  5 


LATITUDE  (DEG) 

(■«.  6 - The  mean  zonal  wind  distribiution  using  the  Rayleigh  friction  profile  given 
hy  equation  ( 14)  and  the  photochemical  acceleration  of  Newtonian  cooling  given  by 
equation  ( 1 2), 


LATITUDE  lOEGl 


(b) 

( If.  7 - The  mean  zonal  wind  diitcibuiiun  computed  at  in  Pig.  S with  the  addition  of 
tropotphertc  forcinf  for  winter  in  the  Northern  Hemuphere.  Part  a.  u part  b. 

f(K);  part  c.  ? (ms"*);  and  part  d,  w * TO  x 10*  cm  (cm  i"*l. 
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latitude  (DEGi 

Id) 

Fif.  - The  mean  zonal  wind  distribution  computed  as  in  Fig.  5 with  the  addition  ol 
tropospheric  forcing  for  winter  in  the  Northern  Hemisphere.  Part  a.  i7  Ims*^):  part  b. 
riKl;  part  c.  Fims'^;  and  part  d.  w X 7.0  < 10*  cmicm  s*') 

3 


LATITUDE  (DEG) 


(a) 


Fig.  8 - Same  as  Figure  8 except  tropospheric  forcing  for  wmmer  in  the 
Northern  hemisphere  is  used.  Part  a,  u tms'*);  part  b w (cm  s 


LATITUDE  (DEGI 
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(b) 


Fig.  10  - The  mean  circulation  of  the  upper  atmosphere  in- 
j eluding  zonal  harmonic,  m • 1;  otherwise  same  as  Figure  8, 
Part  e shows  the  vertical  structure  of  the  wave  amplitude; 
units  are  1 x 10*  cm*  s'*  or  about  10  gpm. 


LATITUDE  (DEG) 


Tiy.  12  - Three  ozone  density  profiles  icm'^)  used  to  compute  the  hcatinii  rate  and 
resultant  static  stability.  ProPile  II  was  used  tot  all  calculations  shown  in  previous 
fifures. 
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Fi|.  13  - Same  as  Figure  ’ but  for  equinox  conditions.  Part  e 
IS  perturbation  heating  rate  in  K day** 
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